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ABSTRACT 

We report the discovery of possible overdensities of galaxies at z ~ 3.7 in Chandra Deep Field 
South (CDF-S). These overdensities are identified from a photometric redshift-selected sample, and 
the BUz-selected sample. One overdensity is identified in the proximity of two active galactic nuclei 
and Lyman break galaxies at z = 3.66 and z — 3.70 at 7a significance level. The other overdensity is 
less significant. It is identified around six z spec ~ 3.6 galaxies at 3er significance level. The line-of-sight 
velocity dispersions of these overdensities are found to be a v ~ 500 — 800 km s -1 , comparable to the 
velocity dispersions of clusters of galaxies today. Through spectral energy distribution fitting, we find 
~ 15 massive galaxies with M > 10 11 Mq around the z ~ 3.7 overdensity. The mass of the z ~ 3.7 
overdensity is found to be a few xl0 14 M Q . Our result suggests that high-redshift overdense regions 
can be found in a supposedly blank field, and that the emergence of massive structures can be traced 
back to redshifts as high as z ~ 3.7. 

Subject headings: cosmology: observations - galaxies: clusters: general - galaxies: evolution - galaxies: 
high-redshift 



1. INTRODUCTION 

In hierarchical galaxy formation, galaxy clusters grow 
through gravitational attraction of matter around high 
a peaks, and they are predicted to be rare at high red- 
shift. Therefore, constraining the number density of 
high-redshift clusters is an important method of testing 
hierarchical models (Evrard et al. 2002; Mantz et al. 
2008). 

Many previous studies have been carried out to 
search for evolved structures at high redshift. Ra- 
dio galaxies and active galactic nuclei (AGNs) are 
considered to be signposts for high-redshift proto- 
clusters, and se veral studies have found proto-clusters 
at 2 < z < 4 (iLe Fevre et al.1 H99l iKurk et all l2000t 



Pentericci et all 120001: iWold et all 120031: iK aiisaw a et all 
20061 : iKodama et al.ll2007l : IVenemans et alJl2007f ). Over- 
dense regions are also identified out to z ~ 6 using galax- 
ies selected with the Lyman break technique or narrow- 
band imaging (e.g, Ouchi et al. 2005). 

When identifying overdense regions, it is necessary to 
compare the overdensity with a control field which lacks 
notable high-redshift overdense regions. For such a pur- 
pose, many studies have used the Chandra Deep Field 
South (CDF-S) due to the wealth of the deep muiti- 
wavelength data over a moderately large field of view 
(|Kaiisawa et al.ll2006t IKodama et al.ll2007t ). However, a 
matter of concern is whether the CDF-S can really be 
considered as a field devoid of high-redshift overdense 
regions. There exist a fair number of high-redshift X-ray- 
detected AGNs in the CDF-S, which may harbor high- 
redshift clusters. 

Motivated by this, we searched for signs of overdensity 
around two AGNs at z = 3.66 and 3.70 separated by 2'. 7 
in the CDF-S, and we report the discovery of possible 
overdense regions at z ~ 3.7 in the CDF-S. 

Throughout this Letter, we assume a cosmology with 
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nonzero cosmological constant (Im et al. 1997), f2 m = 
0.3, fl A = .7 and h = 0.7 cons istent with the WMAP 
cosmology ()Spergel et al.l |2007| ). All magnitudes are 
given in the AB system. 

2. SAMPLE SELECTION OF HIGH REDSHIFT GALAXY 

We used two galaxy samples to identify overdense ar- 
eas at z =3.6-3.7; a -BUz-selected sample, and galaxies 
selected with photometric redshifts. 

2.1. BVz Selection 

The Lyman break technique has proved to be an ef- 
ficient method to select galaxies at high redshift from 
multicolor optical data (Steidel et al. 1996). For exam- 
ple, Giavalisco et al. (2004b) selected B-band dropouts 
at z ~ 4 from GOODS data. However, since the -B435- 
dropout selection contains galaxies over a broad redshift 
range, we defined a color space optimized to select galax- 
ies at z ~3.6-3.7 using available spectroscopic redshifts 
including 17 at z -3.6-3.7 (Vanzella et al. 2006, 2008; 
Popesso et al. 2008). Figure 1 shows our selection 
box, overlaid on galaxies in the CDF-S field. The figure 
demonstrates that our selection criteria effectively filter 
out objects at z < 3.55, and at z > 3.75. For example, 
there are 18 galaxies at z spec ~ 3.47, but only one of them 
makes it into the BV z-selected sample. For the photo- 
metric data, we used the version rl.l Advanced Cam- 
era for Surveys (ACS) multiband catalog of the GOODS 
team (Giavalisco et al. 2004a), and the magnitude cut 
of -Z850 < 26.5 mag was chosen. 

With this method, we selected 245 objects over the 
CDF-S area of —160 arcmin 2 corresponding to a surface 
density of E = 1.53 ± 0.06 arcmin^ 2 (Poisson noise). 

2.2. Photometric redshifts 

The BVz selection is biased for UV-bright galaxies. 
As a complementary sample, we also selected if-band- 
limited [K s < 23.8 mag), z —3.6-3.7 galaxy candidates 
using photometric redshift. 
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We estimated the photometric redshifts of the K s - 
band-limited objects using the Bayesian photometric red- 
shift estimation (BPZ; Benftez 2000). The photometric 
data include 

UU'B4^Veoeir75Zs5oJHK bands, and Spitzer IRAC 
photometry (Giavalisco et al. 2004a). The Near-infrared 
(NIR) photometry was performed using the version 2.0 
released images in the J-, H-, and i^-bands of the ESO- 
GOODS team (J. Retzlaff et al. 2008, in preparation). 
The depth of the VLT/ISAAC data reaches K = 24.7 
at the 5er AB limits, and the simultaneous coverage in 
JHK S bands is 160 arcmin 2 after trimming the edges. To 
create a catalog, we ran SExtractor (Bertin & Arnouts 
1996) in double-image mode, using i^-band detections 
as a reference, and performing photometry in the J- and 
-ff-bands at the K s positions. Auto-magnitudes were 
taken for the J,H , and K s photometry. The IRAC pho- 
tometry catalog was kindly provided by R. Chary (2008, 
private communication). In most cases, we found IRAC 
detections for the -fT-band selected objects. 

To derive photometric redshifts, we used 742 spectral 
energy distribution (SED) templates with a single burst, 
or a constant star- formation rate (SFR), or an exponen- 
tially declining SFR in the form of SFR oc exp(t/r) for 
r = 1, 2, 3, 5, 15, or 30Gyr. These template SEDs were 
generated using a stellar population synthesis model of 
Bruzual & Chariot (2003). The metallicity was varied 
to be OAZq, IZq, or 2Z Q , and we adopted the Salpeter 
initial mass function extending from 0.1 to 100 Mq. The 
SED ages range from 100 Myr to 10 Gyr or to the age of 
the universe at the corresponding redshift. The redden- 
ing parameter Ay was allowed to change from to 1.8 
according to the simple two-component model of Chariot 
& Fall (2000). 

From the comparison of z spec versus z p hot, we find that 
contamination from low-redshift interlopers is minimal 
when selecting z ~3.6-3.7 objects, at the redshift inter- 
val of 3.45 < z p hot < 3.85. When we restrict the redshift 
interval as above, the scatter in the z p hot versus z spec re- 
lation is about 0.1 after excluding several clear outliers. 
Therefore, we choose 3.45 < z p hot < 3.85 as the inter- 
val for the photo- z sample selection. This gives 60 ob- 
jects in total, including 17 spectroscopically confirmed 
galaxies with good z spec quality flags (B or better) at 
z ~3.6-3.7. Note that we use z spec instead of z p hot when 
available, but restrict the spectroscopic redshift range 
to be 3.585 < z < 3.706. With the photo-z selection 
criteria, 73% of z spec =3.6-3.7 galaxies are detected as 
Z =3.6-3.7 galaxy candidates. The contamination rate is 
26%, but the contamination due to low redshift objects 
(z S pec < 3.3) is zero. Most of the contamination comes 
from z spec ~ 3.5 objects. 

3. RESULTS 

We searched for the overdensities by running circu- 
lar top-hat filters with various diameters (0'.5 to 4'. 8, 
corresponding to 0.2-2 Mpc at z = 3.7) on the two- 
dimensional spatial distribution of galaxies in the above 
samples. The number of objects was counted around a 
given point within a circle with a particular diameter, 
producing a density contour map. Figure 2(a) shows the 
surface number density contour map and the location of 
-B^z-selected sample and photo-z selected objects, re- 
spectively. Here, the contours are created using 2'. 4 di- 
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Fig. 1. — V - z vs. B — V colors of objects in the CDF-S 
with z < 26.5 mag. The parallelogram in the middle of the figure 
represents the selection criteria of our SVz-selected sample. This 
is similar to the ZJ-band dropout selection box (thin line; Giavalisco 
et al. 2004b). The spectroscopic data are overplotted. 



TABLE 1 

Significance of overdensity 



sample 


mean density 


1-cr 


A (South) 


B (North) 




arcmin - 2 




deviation a 


deviation a 


By^-selected sample 


1.48 


0.64 


4.94 


2.52 


3.45 < z phot < 3.85 


0.37 


0.27 


6.86 


1.96 


Sum of two samples 


1.68 


0.69 


6.18 


2.37 



ameter filters, which corresponds to ~1.0 Mpc in phys- 
ical size at z = 3.7. The 1 Mpc size is found to be the 
optimal size to identify the z = 3.7 overdensity (see be- 
low). The spectroscopically confirmed galaxies are plot- 
ted with squares (3.45 < z spec < 3.55), diamonds (3.55 < 
Zspec < 3.65), and triangles (3.65 < z spec < 3.75), and 
the open stars represent AGNs at z spec =3.66 and 3.70 
from Szokoly et al. (2004). One density peak is apparent, 
located in the southern part of the CDF-S (A in Figure 
2(a)), around z = 3.66, z = 3.70 AGNs, and (z spec )=3.70 
Lyman break galaxies (LBGs, diamonds). Another pos- 
sible peak exists in the north around (z spec ) — 3.60 LBGs 
(B), but its significance is much weaker. 

To quantify the statistical significance of these peaks, 
we measured the overdensity factor as a number of a 
from the mean surface number density by performing a 
Gaussian fit over the surface number density histogram. 
The significances of the overdensities are summarized in 
Table 1 . The significance of the overdensity is 5 — la for 
the A overdensity, demonstrating that its significance is 
very strong. The significance of the overdensity does not 
change for the top-hat filter size between 0.75 and 1.3 
Mpc. Beyond these sizes, the significance decreases. On 
the other hand, the B overdensity has a significance of 
only 2 — 3(7. 

With the limited number of spectroscopic redshifts, 
assigning redshifts of the overdense regions is not an 
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Fig. 2. — (a, Left): Spatial distribution of the B^z-selected sample (far-left) and photo-z selected objects, both indicated by crosses. 
Squares: 3.45 < z spec < 3.55; diamonds: 3.55 < z spec < 3.65; triangles: 3.65 < z Bpec < 3.75; stars: AGNs at z 3 pec = 3.66 and 3.70. The 
possible overdensities are identified with dashed circles, and their peak locations are at 03 h 32 m 26 s , -27°51 m 06 s (J2000) for the BVz- 
selected sample, and 03 ,l 32 m 24 s , — 27°50 m 50 s , for the photo-z selected objects, (b, Right): Redshift distribution for the spectroscopic 
sample at 3.45 < z < 3.75. The shaded region shows the objects located within 1.0 Mpc radius of the southern overdensity peak, A. 
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Fig. 3. — Massive galaxies (M a tar > 10 11 Mq) in the CDF-S. Left: The optical to 8^tm SEDs for two representative massive galaxies 
associated with the overdense region A. The best-fit Bruzual & Chariot (2003) models to the template SEDs (solid line) are shown. Right: 
the spatial distribution of massive galaxies (square) on the sum of two samples (cross) peaked at 03 h 32 m 23 s , -27°52 m 19 s . 

fields. For this, we examined only the BVz-selected sam- 
ple, since the photo-z sample does not have comparative 
samples readily available. 

We studied the number density of -BVz-selected 
objects in the GOODS-North field to compare with 
the CDF-S. Using the version rl.l ACS multiband 
source catalog, we find 220 i?Uz-selected objects in the 
GOODS-North field (see Giavalisco et al. 2004b). This 
gives the surface density E = 1.38 ± 0.07 arcmin -2 , 
comparable to the number in CDF-S (S = 1.53 ± 0.06 
arcmin -2 ). 

We also looked for a sign of overdensity at z ~ 3.7 us- 
ing the BV z-selected sample in the GOODS-North field. 
We find the strongest overdensity has the overdensity sig- 
nificance of 3.7a. This suggests that the GOODS-North 
may contain a weak overdense area at z ~ 3.7, but its 
significance is much lower than the overdense region we 
identified in the CDF-S. 

4.2. Stellar Population 

In order to find massive galaxies associated with the 
overdense regions, we performed SED fittings to 60 galax- 



easy task. Nevertheless, we find several redshifts to be 
plausible. The redshift distribution of CDF-S galaxies 
with z spec shows three distinct peaks at z spec ~ 3.47, 
Zspec ~ 3.60, and z spec ~ 3.70 (Figure 2(b)). The major- 
ity of z S pec ~ 3.70 objects are located around the over- 
dense area A (filled histogram), as well as two AGNs at 
Zspec = 3.70 and 3.66, suggesting strongly that this over- 
density is at z ~ 3.7. There are also several z spe c ~ 3.47 
galaxies near the A overdensity. However, the overden- 
sity is more likely to be at z ~ 3.7, since the overdense 
area is clearly identified in the contour map of the BVz 
sample which effectively filters out z < 3.55 galaxies. 

Around the weak, overdense region B, we find six LBGs 
with (z) = 3.60, and find a one-dimensional velocity dis- 
persion of a v ~ 500 km s _1 for these objects. The B 
overdense region may be a loosely bound overdensity at 
2 = 3.6. 

4. DISCUSSION 

4.1. Cosmic Variance 

We examined whether the number of high-redshift can- 
didates is abnormally high in this field compared to other 
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ies in the photo- z sample following the procedure of Shim 
et al. (2007), with the SED templates used for the photo- 
metric redshift determination. Note that we have limited 
the SED fitting to objects with IRAC data in order to 
reduce uncertainties and degeneracies arising from using 
only the rest-frame optical data. 

In Figure 3, we show the best-fit SEDs of two repre- 
sentative massive galaxies at the overdense region A in 
the CDF-S. We find 35 galaxies with M star > 10 11 M 
over the whole CDF-S area. More than 43% (15) of these 
massive objects are located within ~ 1 Mpc radius of the 
overdensity peak of the A region. This offers a further 
support for the existence of z ~ 3.7 overdensity. 

We also checked the usefulness of the H — K color selec- 
tion of z ~ 3.7 as advocated by Brammer & van Dokkum 
(2007). They suggest that objects with H-K > 0.9 mag 
are likely to be at (z) — 3.7 based on their photometric 
redshift analysis. However, our analysis shows that only 
about 30% of H — K > 0.9 objects are at z > 3, for both 
the z spec and the z p hot samples. Therefore, we conclude 
that the usefulness of the H — K > 0.9 selection method 
is limited. 

4.3. Overdensity Mass 

We estimate the overdensity mass in several different 
ways for the A area. 

First, we assume that it is a virialized structure, fol- 
lowing Biviano et al. (2006). This assumption is not 
plausible, considering that the age of the universe is only 
1.6 Gyr at z = 3.7, while it takes nearly several Gyr for 
galaxies to cross the 1 Mpc structure at the given ve- 
locity dispersion. Nevertheless, the method can place a 
useful upper limit on the overdensity mass. For the clus- 
ter size, we use the harmonic mean radius of the photo- z 
selected objects within a 1 Mpc radius (twice the filter 
size that gives the maximal overdensity signal) from the 
overdensity center. This gives a size of 0.55 Mpc. As- 
suming that the one-dimensional velocity dispersion, a v , 
of several possible proto-cluster members with z spec is a 
representative value, we estimate the cluster virial mass 
using the gapper estimator (Beers et al. 1990). We find 
that a v — 800 km s _1 by excluding a small redshift peak 
at z spec — 3.66. The derived cluster virial mass is, then 
M mr ~5x 10 14 M o . 

Next, we apply a procedure similar to that of Vene- 
mans et al. (2005) on the z p h ot -selected sample. Suppose 
that the cluster is roughly spherical, and has a 0.55 Mpc 
radius equivalent to the harmonic mean radius as derived 
above. Using the range 3.45 < z < 3.85 as the selec- 
tion window, we find an overdensity factor of S ga i ~ 700. 
Adopting the bias parameter of b — 4 of B-dropouts in 
the CDF-S (Lee et al. 2006), we obtain a cluster mass of 
4 x 10 14 M . 

Finally, we add up the stellar mass of the member 



candidates within the 0.55 Mpc radius, and convert it 
to the total mass assuming a M sta r-^o-Mh a io ratio of 
0.026 - 0.056 of early-type galaxies (Jiang & Kochanck 
2007). Here, the overdensity mass can be expressed as 
M=(N g xfx (Af ator ))/(0.026 - 0.056) M , where N g 
and (M st ar) are, respectively, the number and the av- 
erage mass of galaxies with M star > 10 11 M Q , and / is 
the ratio of the mass of galaxies with M > 10 11 M to 
those with M < 10 11 M . In our case, we have N g = 6 
and (M star ) = 5 x lO n M . We adopt / ~ 4, which is 
the ratio of the luminosity density of L > L* galaxies to 
that of L < L* galaxies in the Schechter function with 
faint-end slope a ~-1.24 (e.g, Rines & Geller 2008). This 
gives an overdensity mass of (2.1-4.6) x 10 14 M 

It is interesting that several independent mass esti- 
mates give the overdensity mass to be M ~ a few 
x 1O 14 M . The implausible assumption of the virializa- 
tion even gives a mass consistent with other estimates. 
This suggests that the virialization of the proto-cluster 
may soon be completed. However, the expected number 
density of halos as massive as a few x 1O 14 M at z ; ~ 3.5 
is extremely low with 10~ 10 Mpc" 3 ()Park fe Kimll2007t 
iReed et aT1l2007t ). and the probability of finding such a 
halo in the CDF-S is only 10" 5 . This structure can be a 
part of the filaments and walls at high redshift as well. 
The existence of another weak overdense area such as the 
B overdensity supports such an idea. A careful analysis 
of simulation data should help us understand the nature 
of these kinds of overdensitics. 

5. CONCLUSIONS 

By examining the two-dimensional distribution of 
3.45 < z < 3.85 candidates in the CDF-S field selected 
from the photometric redshift method or the BVz color- 
color space, we find plausible associations of the over- 
dense region of such objects with AGNs and LBGs at 
Zspec ~ 3.7. The significance of the overdense regions is 
found to be 5 — 7a. The overdense area is abundant with 
massive galaxies, adding further support that this is a 
proto-cluster. The derived mass of the proto-cluster is 
found to be a few x 10 14 M . The existence of the mas- 
sive proto-cluster in the CDF-S shows that one must be 
cautious when examining the significance of high redshift 
overdensities in other fields with respect to the CDF-S. 
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